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Abstract 


Recent work has emphasised the importance of vector and tensor contribu¬ 
tions to the large scale microwave anisotropy fluctuations produced by cosmic 
defects. In this paper we provide a general discussion of these contributions, 
and how their magnitude is constrained by the fundamental assumptions of 
causality, scaling, and statistical isotropy. We discuss an analytic model which 
illustrates and explains how the ratios of isotropic and anisotropic scalar, vec¬ 
tor and tensor microwave anisotropies are determined. This provides a check 
of the results from large scale numerical simulations, conhrming the numerical 
finding that vector and tensor modes provide substantial contributions to the 
large angle anisotropies. This leads to a suppression of the scalar normalisa¬ 
tion and consequently of the Doppler peaks. 


I. INTRODUCTION 

The idea that the breakdown of some fundamental symmetry and the consequent field 
ordering might be responsible for structure formation in the universe is an attractive one. 
Recently we have performed the first complete calculations of the power spectra of perturba¬ 
tions in symmetry breaking theories, including global cosmic strings, monopoles and texture 
i, i- These calculations revealed that vector and tensor modes give a larger contribution 
to the large scale anisotropies than previously suspected, and that their fractional contri¬ 
butions to the total microwave anisotropy power spectrum are comparable for each theory 
considered (see Figure 1). Simultaneous work on local strings [Q has produced compatible 
conclusions. 
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The main implication of the large vector and tensor contribution on large angular scales 
is in reducing the normalisation of the scalar perturbations, which are responsible for the 
Doppler peaks. Once the vector and tensor contributions are properly included, the height 
of the Doppler peaks are low relative to the large angular scale Sachs-Wolfe plateau |ll]. 

The present paper represents an analytical attempt to explain why vector and tensor 
contributions are substantial on large angular scales, using only the most fundamental prop¬ 
erties of the simplest defect theories, namely scaling, causality and statistical isotropy. We 
illustrate the arguments through comparison with the results of numerical computations [|l| . 


II. CAUSALITY AND ANALYTICITY 

As discussed in p, all perturbation power spectra are determined by the unequal time 
correlator (UETC) of the defect source stress energy tensor 0^,^: 

(0/,^(k, r)0pA(-k, t')) = Cf,^,px{k, r, r') (1) 

where r, t' denote conformal time, and k comoving wavenumber. Note that ( 0 ) is real 
because complex conjugation is equivalent to the replacement k ^ —k. The correlators are 
invariant under this replacement because the statistical ensemble is rotation invariant. 

Causality means that the real space correlators of the fluctuating part of must be 
zero for r > r -|- r' 10. Scaling dictates that in the pure matter or radiation eras oc 

px{kT, kr'), where 0o is the symmetry breaking scale and c is a dimensionless 
scaling function. Finally, Qpi, must obey the equations for stress energy conservation with 
respect to the background metric (see next Section). These provide two linear constraints on 
the four scalar components of the source. Any pair determines the other two up to possible 
integration constants. In the matter era the pair 0 and 0"^ 0 provides a convenient choice, 
allowing an analytical integral solution to the linearised Einstein equations. But for work 
including the matter-radiation transition 0 the pair 0oo and 0'^ is better, because it results 
in the correct redshifting away of all components of the source stress energy inside the 
horizon. In this paper we shall use both pairs - 0 and 0“^ in our analytical discussion of an 
‘incoherent’ model, and 0oo and 0"^ for a numericaly solved ‘coherent’ model. In the former 
case, we shall constrain 0 and 0“^ so that on subhorizon scales the source is negligible (see 
Section IV). 

The unequal time correlator in k space is the Fourier transform of the real space corre¬ 
lator: (0jj(k, r)0fci(—k, r')) = / (i^re“**^'’'(0jj(r, r)0M(O, r')). The integral is hnite because 
the real space correlator has compact support, and it follows that the unequal time corre¬ 
lators are analytic in k for all hnite k. They may thus be expanded as Taylor series in the 
Cartesian components ki about fcj = 0. As ki tends to zero, isotropy and symmetry impose 

him k ^ 0 (0p(k, r)0M(-k, r')) = A6ij6ki + B{Sik6ji + dudjk) (2) 

with A and B independent of k. 

The trace scalar, anisotropic scalar, vector and tensor components of a tensor Tij are 
given by 

T.,(k) = + (kkj - + (kjf + kjT^) + TT 
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FIG. 1. The contributions to the total anisotropy power spectrum from scalar, vector and tensor 
components, in the theories of global strings, monopoles, texture and nontopological texture (taken 
from ref. [1]). 
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FIG. 2. The importance of the long wavelength modes in the anisotropy power spectra from cos¬ 
mic textures. The power spectra due to scalar (dotted), vector (dashed) and tensor (long-dashed) 
components of the sources are compared to those where the source stress energy components ©oo 
and 0'^ as well as the vector and tensor stresses are set zero for all fcr > 5. The upper curves show 
the full spectra, the lower ones the results where the cutoff is imposed. 
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( 3 ) 


where ki = kijk. Expressing the trace T, T'^, and in terms of Tij (see e.g. [Q) one 
hnds that the only nonzero correlators consistent with statistical isotropy and homogeneity 
are (TT), (TT*^), (T'^T'^), (T^T^) and {T^T^i). From (^ one can compute the small k 
power spectra of the anisotropic scalar, vector and tensor stresses. One hnds they are in the 
ratios 


(|0^P):(|0rP):(|0?;f) = 3:2:4 (4) 

where all indices are summed. Thus in a causal theory anisotropic scalar, vector and tensor 
stresses have white noise components at small k with related amplitudes. A similar argument 
shows that the correlator (0oo0p') ~ at small fc, implying that (0oo0'^) vanishes like 
at small k. Likewise (00‘^) vanishes like k"^ at small k. So for either of the two choices 
discussed above, the two scalar source components are uncorrelated outside the horizon. 


III. SUPERHORIZON MODES 

In the cosmic defect theories, perturbations are predominantly produced on the horizon 
scale. Studies show that the unequal time correlators take the predicted white noise form 
for fcr < 5 or so, and decline strongly at larger kr. To the extent that the horizon scale 
modes reflect the causality constraints discussed above, the latter translate into dehnite 
relations between the scalar, vector and tensor perturbation power spectra. In Figures 2 
and 3 we show the CMB anisotropy power spectra calculated in the cosmic global string and 
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FIG. 3. As in Figure 2 but for global strings. 


texture theories respectively, with and without a cutoff where we switch off the source stress 
tensor for kr > 5. The Figures show that in the texture theory the effect of suppressing 
the source for fcr > 5 is relatively minor. For strings, there is a larger effect, but even here 
the ratios of scalar to vector to tensor anisotropies are not much affected. We conclude 
that the contributions from kr < 5, which we shall term superhorizon modes, are certainly 
important in both theories and give at the least a rough measure of the importance of the 
scalar, vector and tensor contributions to the large angle anisotropies. 


IV. INTEGRAL CONSTRAINTS 

The fact that a cutoff on subhorizon scales does not greatly affect the large angle Ci 
spectrum has important implications. It means that the short distance structure of the 
individual defects is not important in determining the qualitative character of the large 
angle anisotropies, such as the relative scalar/vector/tensor contributions. 

Consider the effect of modelling the sources using a ‘smoothed’ ©oo tensor, one where 
we impose a cutoff at kr ~ 5. We feed in the smoothed ©oo and ©'^ into the stress energy 
conervation equations: 

©00 + “(©00 + ©) = !!; 11 + 2—11 = —;y(© + 2©'^). (5) 

a a 6 

where 11 = cIj©oi and © and ©‘^ are dehned in the previous Section. It is straightforward to 
see that the solutions for 11 and © are well dehned. 

In this scheme, we deduce important relations between © and ©'^. Equations (|^) are 
easily integrated to obtain ©oo in terms of © and ©‘^: exchanging the order of the double 
integral we get 

000 = T-"1 r dT'2T'e + - i)(e + 2 e®)(T')|, ( 6 ) 

Jo 6 T' T 
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where we used a(r) oc in the matter era. But as argued, the smoothed ©oo is identically 
zero inside the horizon. It follows that both the and t~^ coefficients integrate to zero. 
The former gives 

roo 1 

/ dr'[2r'0 + -fcV3(0 + 20^)(fc,r')] = 0, (7) 

Jo 6 

and the latter gives 

roo 

/ dT'T^(Q + 2Q^)(k,T') = 0. (8) 

Jo 

This imposes a negative correlation between 0 and , and guarantees that 11 vanishes 
faster than inside the horizon. The constraints and (@) will turn out to be remarkably 
powerful when building models for the superhorizon components of 0 and 0'^. 


V. PERTURBATIONS IN THE MATTER ERA 


We wish to compute the large angular scale anisotropies produced in the matter era. For 
this purpose we use the following integral solution to the linearised Einstein equations in a 
matter dominated universe |^: 


6T 
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h-r = -16vrG^, 




drh,,,o(r, n(ro - 
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tensor 


+ 20'')(r',k) - kyyfQ^r'y)] 

or T 


r vector _ 

'\j,0 - 


Y. To = 16 >rG £ dr'i^Y (,/) 


h^ST’iT,^) = 16t,G £ dT'kh'*lG,{T')G2{T) - G2(T')Gl(T)]e5(T'x) 

cos(fcr) sin(A;r) cos(fcr) sin(fcr) 

= /, xo- r, — ^ 2 {t) = 


(kr) 


{ky ’ 


(fcr) 


(fcr)2 ’ 


( 9 ) 


where Gi and G 2 are the two homogenous solutions to the tensor (gravity wave) equation. 

The model we shall consider is one in which the components of 0*^ have the following 
unequal time autocorrelators: 


(167rG)^(0(k, r)0(—k, t')) = 6{e — kT)5{T — t')A 
(167rG)^(0‘^’^’'^(k, r)0(—k, = 6{e — fcr)5(r — (10) 

where 6 is the Heaviside function, and we dehne (0^(r)0^(r')) = |(0^(T)0h(r')), and 
(0'^(r)0^(r')) = |(0f {t)QY all indices summed. 

The sources are nonzero only on ‘superhorizon’ scales {kr < e) and they are uncorrelated 
except at equal times. This latter property means that the model is ‘totally incoherent’, 
in the terminology of ref. [Q. These correlators are not strictly causal - in real space they 
take the form r“^[sina; — xcosx] where x = er jr - but they are small and oscillatory beyond 
r r-u T for e = 5. So the violations of causality are small. Rotational invariance forbids any 
cross-correlation between scalar, vector or tensor modes. There is however one more allowed 
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cross correlator, namely that between the isotropic and anisotropic stresses. The argument 
given in Section III implies that (00'^) vanishes as for small fc, but it cannot be zero 
because of the constraint (H). We choose to model it as 

(167rG')^(0(k, r)0‘^(—k, r')) = 6{e — A;r)(5(r — r')(—(11) 
If we now compute the equal time correlator of the the constraint (§), we determine 


do 


Similarly we compute the equal time correlator of 


and obtain 
4 




( 12 ) 


(13) 


These equations yield = —2.47^'^ and A = 4. 07 for e = 5. At kr = 5, there is 
a mild inconsistency with the bound < VAA^^ so we shall adopt = —2A^ and 

A = 4^5. 


VI. THE DELTA FUNCTION APPROXIMATION 


The procedure is simple in principle: the correlators (|10|) translate into correlators of 
the metric perturbations and thus into correlators of the temperature perturbations, equiv¬ 
alent to the anisotropy power spectrum Q. But in order to compute the relevant integrals 
analytically, we shall make two approximations. The hrst is that we shall replace metric 
perturbation unequal time correlators under r integrals using the following formula: 


(A(t)B(t')} (5(r - 


(14) 


The weighting function is chosen so that the integrals Jq^ dr Jq^ dr' of both sides are guar¬ 
anteed to be equal for all Tf. The formula is also invariant under changing variables from 
r to any other function /(r). The second approximation is to use the fact that the Greens 
functions in (^ fall off strongly with r. This means that the metric perturbations fall off 
rapidly beyond kr = e, which justihes us simply setting them zero beyond that point. 


VII. CMB ANISOTROPIES 

In the usual way we expand the microwave sky temperature in spherical harmonics 
6T/T = I] 0mh)m(6', 0), and compute Ci = (|a/mP)- The formulae for the contribution to 
the integrated Sachs Wolfe effect from trace scalar and anisotropic scalar contributions are 

i: 
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where At = tq — r, tq is the conformal time today and as above, () denotes ensemble 
averaging. The scalar metric perturbation components are given from (^: 


ill = —WttG j dr'(r7r)®(0 + 20'^)(r') 

/i 2 = —IbvrG J dr'{ t'/{t') (16) 


with k dependence implicit. 

The vector and tensor contributions to Ci are 


cr 
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(17) 


where {h^{T)h^{ t')) = i(h^(r)h^(r')), and (h'^(r)h^(r')) = \{hY{T)hY{ t')) , with all i 


in¬ 


dices summed. 

We now compute the relevant metric perturbation correlators: from (pT)), (^) and (pTl), 
using = —2A^ and A = ‘iA^ as discussed above, we obtain 


(fc2(r)=) = 

{Ai(r)fe 2 (r)) « 

{h''(T?) = TtM'' 

{h^{TY) = J dT'G{T,T')‘^A^ 

G{t, t') = fc^r''^[Gi(r')G2(r) - G2(r')G'i(r)] (18) 


where we have evaluated the scalar correlators at fcr = 5, and the tensor modes Gi and G 2 
are given in equation (^). The tensor integral is straightforwardly performed, yielding 

AT 1 nc; 

(h^(r)2) = ^^( (1 _ z^)sm{2z) - 105zcos(2^) + 6/ - 14^^ - 35^^) (19) 

where z = kr. This function is ~ ^^r^/81 at small fcr, identical to the vector expression. 
But for larger kr it is suppressed, with the suppression factor being 0.29 at kr = 5. The 
suppression is due to the oscillatory nature of the tensor compared to the vector response. 

We now compute the integrals in ([T7|), starting with the tensor contribution Gf. The 
delta function allows one of the r integrals to be performed. Then we change variables 
from k to X = kAr. The Heaviside function gives the upper limit x < e(ro/r — 1), or 
T < tq/{ 1 + x/e). Exchanging orders of the integrals we hnd 
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^(/-2)! Jo 


dx n, , fTo/(l+x/e) 


dr d 
{tq - r)3 dr 




( 20 ) 
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FIG. 4. Tensor anisotropy power spectrum as computed in the analytical model presented here. 



For large I the integral is dominated by large x, since ji{x) ~ at small x. Bnt at large x, 
r << To and the r integral is trivial. Thus one finds at large I 


C 


81 271 Jo x^ 


Jii^) 




- 3 \ 


where we have used 


/ 


dxjf{x)x'^ 


TT r(i-n)r(/ + i + f) 
22-r(i-f)2r(/ + |-f)' 


( 21 ) 


( 22 ) 


We have computed the integral for the tensor contribution (^0[) at low I using Mathematica, 
to check that the model reproduces the shape of Z(/ + l)Ci seen in the plots of Figure 1. 
Figure 4 confirms that this is indeed the case. 

The vector integral is performed similarly, to obtain 


Cj' = ^ dxx^[{ii{x)/x)’f[^ - ^ + ln(l + e/x)]. (23) 

where prime denotes differentiation with respect to x. Approximating the square bracket 
with its leading large x behaviour, integrating by parts and using Bessel’s equation for ji(x), 
we obtain at large I 


Cl 


V 




(24) 


The scalar contribution is evaluated using ([T^) and (^), again making the large x approxi¬ 
mation, giving 


^scalar 


= A^- X^dxi—jf + -J- + 


277 


120 


81^ 


217-’ 


2x3 


(25) 
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After integrations by parts and using Bessel’s equations we get 


^scalar 


1 

1248 




(26) 


Now the amplitudes , AA and A^ are related via in the ratio 3: 1: 1. Thus in the 
various approximations we have made, the ratio of the scalar to vector to tensor contributions 
to the large angle anisotropies are 


(.scalar . (jV . (jT ^ . ( . g ^g 


(27) 


which is our main result. The calculation demonstrates the relative importance of the vector 
and tensor modes, consistent with the numerical results shown in Figure 1. Given the crude 
nature of the model used, the agreement is actually surprisingly good. The weakest point 
in the model is that it involves a free parameter e, and the Q’s obtained are proportional to 
e^. It seems plausible that e should be the same for the scalar, vector and tensor stresses, 
but we have not found any argument as to why this should necessarily be true. 

Let us summarise the approximations and assumptions implicit in (^71): 

• We assumed that superhorizon modes with kr < 5 dominate. 

• We modelled the unequal time correlators as delta functions with a horizon scale cutoff. 

• We made the approximation of pure matter domination for the background spacetime. In 
this approximation the Ci spectra obtained are scale invariant at large /. 

• We replaced certain functions with amplitudes times delta functions in order to perform 
the relevant integrals. 

With all of these caveats, we feel that the model provides useful insight into the relative 
importance of scalar, vector and tensor contributions to the large angle anisotropies. The 
model explains why vector perturbations dominate over tensors, and why the combined 
vector and tensor contribution is comparable to that from scalars. 


VIII. NUMERICAL SOLUTION OF A COHERENT MODEL 

As a further model we have considered the case of a completely coherent source in which 
the unequal time correlators of ©oo, 0'^, 0^, and 0'^ are all proportional to the product of 
Heaviside functions, for example setting 

(0oo(k, r)0oo(-k, r')) = (rr')“50(e - /cr)0(e - kr). (28) 

Note that we need to model the white noise contribution, but as mentioned above, the 
cross correlator must vanish at small k. So in this model we will assume that < ©oo©"^ > is 
identically zero, and therefore solve for the ©oo and ©‘^ contributions separately. Note that 
with this choice of variables, the constraints (^ and (@) are automatically satished. 

We have used this model in the full Boltzmann code developed in as usual with e = 5. 
The results are shown in Figure 5. The anisotropic scalar, vector and tensor Q’s have been 
scaled so that the stress tensor white noise superhorizon amplitudes are in the correct ratios 
(i)- kas been scaled with Gf. It is a remarkably small contribution. This may be 
understood by solving the equations for stress energy conservation for 0, whence one hnds 
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FIG. 5. Angular power spectrum of anisotropies generated by a simple coherent model of 
scaling sources, correctly incorporating the superhorizon constraints on the relative importance of 
the anisotropic stress 0'^, vector and tensor 0^ perturbations to the source stress tensor. 

that 0 is actually very small in the coherent model at horizon crossing fcr ~ 5, so that from 
(y) its contribution to the anisotropy is small. 

The coherent model provides a useful comparison to the previous incoherent model. The 
broad agreement between the two models suggests that our main result (^) is actually 
insensitive to the detailed nature of the source. An advantage of the coherent model is that 
we can more easily incorporate the matter-radiation transition, giving rise to departures 
from scale invariance in the Ci spectrum qualitatively similar in character to those observed 
in the realistic source calculations. And as seen in Figure 5 the model gives a reasonable 
impression of the main features of the realistic calculations in Figure 1, at least on large 
angular scales. 


IX. CONCLUSIONS 

In this paper we have developed a set of physically reasonable models for the perturba¬ 
tions generated on superhorizon scales by causal sources. We gave some rigorous and some 
approximate arguments that the large angular scale anisotropies due to the scalar and vector 
plus tensor modes are in general similar in magnitude. If the vector and tensor contributions 
to the large angle anisotropies are large, the scalar normalisation is lower and the Doppler 
peaks due to scalar perturbations are small compared to the large angle Sachs-Wolfe plateau. 

Let us close by mentioning some loopholes in the above arguments, which make it possible 
to circumvent the conclusion that causal sources are unlikely to have large Doppler peaks. 

• If subhorizon modes with kr > 5 dominate the anisotropies then our arguments do not 
apply. Sources in this category have been explored by Durrer and Sakellariadou |P]. 

• If the shape of the unequal time correlators, parametrised in our models by e is different 
for the scalar, vector and tensor components, then the Ci contributions could be strongly 
affected, since Ci oc e^. One could imagine a model where e for scalars was larger than for 
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vectors and tensors, bnt even here one would probably not find sharp Doppler peaks, since 
increasing e is likely to increase the incoherence of the source and thus smoothe out the 
Doppler peaks (this is apparent in Figures 2 and 3). 

• One could consider sources like those in [Q in which the anisotropic stresses are by 
construction zero outside the horizon. In such a model the superhorizon constraint (^) is 
satisfied with all terms being zero. In the model of this is true because the real space 
stress energy master functions were taken to be spherically symmetric, clearly a special case. 

• We have assumed perfect scaling of the sources, and matter domination. There is some 
violation of scaling due to the matter-radiation transition, but this is a small effect on large 
angular scales. Stronger departures from scaling would result from a non-minimal-coupling 
mass term R9^ for the Goldstone bosons. We are currently exploring this possibility. 
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